Dynamical scaling and intermittency in shell models of turbulence 
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We introduce a model for the turbulent energy cascade aimed at studying the effect of dynamical 
scaling on intermittency. In particular, we show that by slowing down the energy transfer mechanism 
for fixed energy flux, intermittency decreases and eventually disappears. This result supports the 
conjecture that intermittency can be observed only if energy is flowing towards faster and faster 
scales of motion. 



Turbulent flows show very different behaviour at 
changing the embedding physical dimension. The most 
spectacular change is in the reversal of the energy flux, 
from large to small scales in 3d and viceversa in 2(i 
Moreover, the three-dimensional forward energy flux is 
strongly intermittent while the two-dimensional inverse 
energy transfer is almost Gaussian The presence of 
strong anomalous fluctuations in the 3d flux is believed 
to be connected to the existence of a hierarchical or- 
ganization in the eddy turn over times, r^, at different 
scales. Simple dimensional arguments predict the eddy 
turn over time of velocity fluctuations at scale r in the 
inertial range to be of the order of ~ r^/^ . For three 
dimensional turbulence, the scenario is the one of an en- 
ergy cascade from slow (large eddies) to fast (small ed- 
dies), without the possibility for fluctuations at different 
scales to equilibrate. This mechanism should be at the 
origin of the burst-like structure in the forward energy 
transfer mechanism: any unusual fluctuation in the en- 
ergy content at a given scale propagates to smaller and 
smaller scales until it is re-adsorbed by viscous mecha- 
nism. Hence, this is the cause of strong inhomogeneity 
in the spatio-temporal statistical properties of the en- 
ergy transfer process. The same qualitative arguments 
capture the absence of intermittency in the inverse 2d 
process, where energy flows from fast to slow modes, al- 
lowing the receiving modes to feel only the mean fluc- 
tuations of the unstable "mother" eddy. Intermittency 
in the related problem of scalar passive/active quantities 
advected by turbulent flows has been connected to the 
existence of a "dissipative anomaly" , i.e. the presence 
of a non- vanishing scalar dissipation in the limit of small 
molecular diffusivity. Intermittency in the latter case has 
a natural Lagrangian interpretation connected to the ex- 
istence of particles which separate even if starting in co- 
inciding points 0, 3. Such a mechanism is absent in 
inverse cascade, pointing toward the conclusion that in- 
verse cascade regimes cannot be intermittent. Here we 
intend to investigate a similar issue on a purely Eulerian 
base. 

Intermittent fluctuations are usually related to the 
scaling properties of the Navier Stokes equations . Let 
us denote hy v{x, t) the velocity fleld of a turbulent flow 
satisfying the Navier Stokes equations: 

dtv + V •Wv = -Vtt + lyAv + f (1) 



where tt = p/p {p being the pressure and p the den- 
sity), v is the kinematic viscosity of the flow and / is 
the (large scale) external forcing. Equations ^ show 
the remarkable property to be invariant under the scale 
transformation 

Xx , X'^v , X^-^t , V A^+'V. (2) 

Note that ((2Jl implies for the energy dissipation: e 
^3h-ig^ i.e. e is constant if ft, = 1/3 as predicted by the 
K41 theory. In the multifractal theory of turbulent flows, 
h is supposed to be a fluctuating quantity, although on 
the average the rate of energy dissipation is constant Q . 
The above scenario is complemented by the assumption 
that the statistical fluctuations are described by a scaling 
invariant probability distribution ^ X^'^^^'^^ where 
D{h) can be interpreted as the fractal dimension related 
to the fluctuations A''. Regardless of the geometrical in- 
terpretation of D{h), the basic physical question is which 
is the mechanism determining the fluctuations of h and 
why there are fluctuations in the energy flux. The above, 
scale invariant, scenario of the non equilibrium statistical 
properties of turbulent flows suggests a simple although 
not trivial picture of intermittency, related to the energy 
cascade mechanism. At scale r, the amount of kinetic 
energy due to turbulent fluctuations can be estimated as 
{SrvY, where 5rV = v{x + r) — v{x), and where we neglect 
vectorial indexes for simplicity. We expect that the rate 
of energy flux at scale r, denoted by e^, is of the order of 
{6rvY /Tr- Because energy transfer is due to non linear 
interactions, Tr can be estimated as r / 5rV. Thus we ob- 
tain {Srv)'^ ^ re(r) which, as expected, is invariant under 
(0). As a result, the energy transfer statistics is strongly 
correlated to the fluctuations on the energy contents at 
different scales, {6rv)'^, and to their dynamical proper- 
ties, Tr- Moreover, the presence of strong intermittent 
fluctuations in the energy contents at different scales, re- 
flects, via the equation of motions, into non-trivial fluctu- 
ations of the local eddy turn over times. Indeed, previous 
theoretical and numerical works have demonstrated that 
spatial and temporal properties of the energy cascade 
mechanism are strongly correlated 0, . 

The previous phenomenological arguments suggest 
that intermittency can be observed only if energy is 
flowing towards faster scales of motion. It is therefore 
tempting to argue that by decreasing the scaling ex- 
ponents of the eddy turnover times along the cascade, 
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i.e. by slowing down the energy transfer mechanism 
for fixed energy flux, intermittency should decrease and 
eventually disappear. Our aim in this letter is to provide 
clear evidence that the above conjecture holds. 

Our analysis will be performed in the framework of 
shell models of turbulence (see P, Q, 13 and references 
therein). The motivation to use shell models as a possi- 
ble surrogate of the Navier-Stokes dynamics is twofold. 
First, shell models proved to be very successful in re- 
producing many of the statistical features of both 2 and 
3 dimensional turbulent flows, being at the same time 
much easier to simulate numerically. Second, they are 
flexible enough to allow a structural change in their equa- 
tions of motion which will allow us to directly probe the 
importance of time dynamics in fixing the intermittent 
properties of the energy transfer process (see below) . 

In a shell model, the basic variable describing the 've- 
locity field' at scale r„ — A~"ro = fc~^, is a complex 
number u„ satisfying a suitable set of non linear equa- 
tions. There are many version of shell models which have 
been introduced in literature (see for a recent review) . 
Here we choose the one proposed in which is an im- 
proved version of the so-called GOY model 0, ITU 

-vk\un -t- /„ (3) 

where A = 2, c = — (1 -I- 6) and /„ is an external forcing. 
In shell models, we can associate E^v to u„. Clearly equa- 
tion (O satisfies the scaling (|2Jl. The important point on 
shell models like Q is that the statistical properties of 
intermittent fiuctuations, computed either using u„ or 
the instantaneous rate of energy dissipation, are in close 
qualitative and quantitative agreement with those mea- 
sured in laboratory experiments, for homogeneous and 
isotropic turbulence. Thus, shell models provide an use- 
ful tool to investigate in a simple way the physical con- 
sequences of scaling ^ and intermittency. Moreover, at 
variance from 3c? Navier-Stokes equations, the compu- 
tational complexity grows with Reynolds only as i?e^/^, 
allowing for reliable numerical studies also at very high 
Reynolds numbers. 

It is easy to realize that also in shell models, T(fc„) goes 
as as predicted by i.e. the characteristic time 

for energy transfer decreases quite fast as the scale is 
decreasing. If t(A;„) is independent of the constant- 
flux argument would suggest energy equipartition among 
all scales, i.e. (Srv)"^ ~ e^- To be more specific, by con- 
straining ourselves on shell models, let us modify © in 
such a way that there still exists an average Reynolds in- 
dependent rate of energy dissipation e — {i^J2n ^n\'^n\'^) 
while the statistical properties are invariant under the 
scaling transformation: 

kn X'^k„ , Un A''u„ , t X'^'H , 1^ X^+'^v 

(4) 

where cc is a fixed parameter in the equations. Note that 
Q to hold, we should require that the non linear terms 



scale as A^^"^. Note also that e ~ iy{Vv)'^ ~ X^''-'^, thus 
we do not change the constrain on the energy flux. 

In order to satisfy the scahng we introduce the 
following equations: 

dUn . , i^x 1 xr * I J, * 

— = ikg A:,Ju„+iU„+2 + 6u„_iU„+i 

+ ^-y^W,i-2Wn-l] - vkl^'^kl^^'Un + fn (5) 

where ko is the largest scale in the system. Clearly, we 
obtain the old shell model for x = 1, while a direct inspec- 
tion of © shows that (01 is satisfied for any x. Moreover, 
for 1^ = 0, the generalized energy, Q = fcg~^E„A:^~^|M„p, 
is conserved by non liner interactions. Thus, in a statis- 
tically stationary regime we have the 'energy' budget: 

n 

where P oc Real{{J2n f^n^^'^nfn)) is the energy input. 
Note that in this shell model energy fluctuations at scale 
kn are of the order Q„ — A:q~^A;^~^|m„P while the eddy- 
turn-ovcr time for the energy transfer is of the order 

T(fc„) ^ l/(fc^ifc>„). (6) 

Thus the energy flux through wavenumber kn can be 
proved rigorously to be of the order of n„ = A:„|u„p 
as for the original shell model. Therefore, for all x, by 
keeping (n„) = const, we can achieve a constant energy 
flux: 

{\un\')-k-\ (7) 

In summary, while for all x the shell model (O produces 
the same average energy flux, its dynamics is character- 
ized by a different time-scaling properties which changes 
the eddy-turn-over hierarchy in the inertial range. In par- 
ticular for small x, say x ~ 1/3, by combining @ and 
{Tj) we should expect all eddy turn over times to be of 
the same order in the inertial range and therefore the en- 
ergy to reach a quasi-equipartion state for flxed value of 
e. The shell model provides us the interesting possibility 
to study intermittency as function of x, i.e. as a function 
of the dynamical scaling I^J and © . Intuitively, one can 
imagine that decreasing x from its Navier-Stokes value 
X = I induces a smoother and smoother energy transfer 
process towards smaller scales, fluctuations among dif- 
ferent scales tend to equilibrate each other and, conse- 
quently, non-Gaussian fluctuations are depleted. Thus, 
as a function of x the shell model |(SJ) should exhibit a 
kind of phase transition from 'strong' intermittent fluctu- 
ations at X = 1 to Gaussian non intermittent fluctuations 
at X small. In this paper we support the above conjecture 
by numerical simulations of the model (O. Moreover, we 
will give a theoretical argument to estimate the value of 
Xc below which intermittency starts to be depleted. 

In order to keep e constant, regardless of the value of 
v and x, we use the forcing /„ = A„/u* for n = 1,2 



and /„ = for n > 2. We have performed a se- 
ries of numerical simulations for different values of x, 
keeping the parameter b = —0.4 fixed. As a measure 
of intermittency, we compute the generalized kurtosis 
G2p{kn) ~ (|u„|^P)/(|u„|^)P. Figure (P) shows the value 
of 6*4 (fcn) for different values of fc„ as a function of x, 
while in the inset we show the value of Ge{kn)- As one 
clearly see, for x smaller than 0.5 a sharp decrease of G4 
and Gq is observed. For small values of x both G4 and 
Gg become equal to their Gaussian value. The above 




FIG. 1: The value of G2p(fc„) = iM'^") / {[unf}" for p = 2 
and p = 3 (in the inset) as function of x and for two different 
values of kn, namely n — 10 and n — 15. The predicted value 
for the Gaussian Statistics (horizontal dotted line) is reported 
for comparison. 

results tell us that intermittency is depleted for small x 
in agreement with our intuitive argument. Also, within 
the numerical error bars, for x > Xc ^ 0.5, intermittent 
fluctuations seem to weakly depend on x. 

We want now to understand why the transition to no- 
intermittency fluctuations is observed for relatively large 
value of X, namely for a; ~ 0.4. As we shall see in the fol- 
lowing, Xc — 1/3 + Ax where Ax is due to intermittency 
correction. The argument goes as follows. Let us consider 
two scales kn and kn+m ■ The corresponding times for the 

energy transfer are T(fc„) - fc^ and T(fc„+™) - fc,'j+„i, 
where h and h are the values of the 'local' fluctuations of 
Un ^ k~^ and Un+m ~ fcn+m respectively. For x = 1 the 
probability that the ratio Xm = T{kn+m)/T{kn) is larger 
than 1 is extremely small, already for m = 1,2. When x 
becomes smaller than 1 the probability P{xm > 1) start 
growing, even for small m. If Xm is much larger than 
1, then the energy transfer from scale fc„ to scale kn+m 
is stopped and energy tends to equipartition. In order 
to estimate Xc, let us compare, for a given scale sepa- 
ration km, the ratio between two eddy turn over times 
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FIG. 2: Probability distribution of log{T{ki4) /T{ki3)) for two 
different values of x, namely x = 1 and x = 0.5. The theory 
presented in the paper predicts that intermittency is depleted 
at a; = 0.5 because T(fc„+i) ~ T(fc„). In the inset we show the 
probability distribution of log{T{ki7) / rikis,)) 

The above expression tells us that Tm{q) > 1 for 
X < Xc{q) = —C{—q)/q- Let us note that, because of 
convexity properties of C('/); ^ciq) is an increasing func- 
tion of q. Thus, there is not a single value of x be- 
low which intermittency is depleted, rather the transi- 
tion to equipartition is continuous, different moments 
of the eddy turn over ratios behave in slightly differ- 
ent way. As a simple guide line the transition region 
is [xc{l),Xc{2)] — [0.4,0.42], estimated by using the 
D{h) curve which fits the Ciq) exponents as given by 
the She-Leveque formula |l2l |. A direct numerical in- 
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FIG. 3: Same quantities as in figure Q for the model with 
parameter b — —0.8. The increase of intermittency at a; = 
1 increases the critical value of x for which thermodynamic 
equilibrium is attained. 



vestigation gives further support to our previous argu- 
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ment. In figure ^ we show the probabihty density func- 
tion of the ratio ■7"(fci4)/T(A:i3), T{kn) being the instan- 
taneous eddy turnover time at shell A;„. For a; = 1, 
we expect that riku) < r(fci3) with probability close 
to 1, while, according to our estimate, we expect that 
for X — 0.5 T(fci4) T^kis) and intermittency is de- 
pleted. This is clearly shown in figure where in the 
insert we also show the probability density function of 
log{T{kn)/T{kiz))- Hence, we can conclude that for all 
scales fc„ the eddy turnover times become almost equal 
producing a kind of thermodynamic equilibrium in the 
generalized energy. 

As a further check of our argument, we consider the 
shell model (jSl for h = —0.8. In this case, the model at 
X = 1 shows larger intermittent correction with respect 
to the case b = —0.4 previously considered. According 



to our argument for the transition to occur, we should 
expect that for b = —0.8, depletion of intermittency takes 
place for larger values of x, which is indeed the case as 
shown in figure Q. 

In summary, we have presented the following main re- 
sults: (i) we introduced a new version of the shell model 
which satisfies a generalization of the dynamical scaling 
(ii) we have proposed a simple, although non trivial, 
argument for understanding how intermittency can de- 
pend on the scaling properties of the eddy turnover time; 
(iii) we have shown, by numerical simulations, that our 
argument is correct; (iv) we have provided a multifrac- 
tal estimate of the critical value, Xc where intermittency 
should disappear. 

We thank F. Toschi for discussion in a early stage of this 
work. 
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